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Abstract 



We describe a web-based tool ( PASCal; principal axis strain calculator, http: / / pascal.chem.ox.ac.uk) 



designed to simplify the determination of principal coefficients of thermal expansion 
and compressibilities from variable-temperature and variable-pressure lattice param- 

(N 

eter data. In a series of three case studies, we use PASCal to re- analyse previously- 

published lattice parameter data and show that additional scientific insight is obtain- 
ed 

able in each case. First, the two-dimensional metal-organic framework Cu-SIP-3 is 
found to exhibit the strongest area-negative thermal expansion (NTE) effect yet 
observed; second, the widely-used explosive HMX exhibits much stronger mechanical 

•i-H 

anisotropy than had previously been anticipated, including uniaxial NTE driven by 
thermal changes in molecular conformation; and, third, the high-pressure form of the 
mineral malayaite is shown to exhibit a strong negative linear compressibility (NLC) 
effect that arises from correlated tilting of SnOg and SiC>4 coordination polyhedra. 
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1. Introduction 

The mechanisms by which crystalline materials respond to changes in temperature 
and pressure are as important for the valuable insight they provide into the funda- 
mental chemistry of the solid state as they are for a range of practical applications. 
Topical examples include negative thermal expansion (NTE) in metal-organic frame- 
works (MOFs) (Wu et al, 2008; Lock et al, 2010), negative linear compressibility 
(NLC) in molecular co-crystals (Fortes et al, 2011) and extreme gating responses to 
guest sorption in microporous materials (Serre et al, 2002; Mulfort et al, 2010; Wal- 
ton et al, 2011). From an experimental viewpoint, crystallographic techniques offer a 
direct and straightforward measure of any structural variation: the evolution of unit 
cell parameters describes the bulk mechanical response of the material, with the cor- 
responding changes in atomic positions indicating the microscopic mechanism respon- 
sible for this response. So, for any given system, determination of the temperature- 
and/or pressure-dependence of its unit cell parameters is a simple but effective means 
of diagnosing anomalous (and useful) mechanical behaviour. 

Because the unit cell is a construct of convention rather than uniqueness, the ques- 
tion as to whether variations in unit cell geometry actually reflect the fundamental 
mechanical response 1 of a material is an important one. Put differently: can we be 
certain that an alternate choice of unit cell — however unconventional — would yield the 
same lattice expansivities and/or compressibilities? The answer, of course, is that the 
values determined nearly always depend on the particular unit cell chosen. A canoni- 
cal example is given by a simple two-dimensional molecular framework [Fig. 1(a)] (a 

lower dimensional analogue, perhaps, of MIL-53 (Serre et al, 2002) or Ag3[Co(CN)e] 

1 In this manuscript we use the term 'mechanical response' to refer to those mechanical properties most 
readily measured crystallographically: hydrostatic compressibilities, the bulk modulus, the pressure 
derivative of the bulk modulus, and the thermal expansivities. In other contexts, the same term is 
sometimes equated with the full elastic tensor and includes e.g. the elastic compliances, Poisson's 
ratio and shear moduli (Nye, 1957). 
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(Goodwin et al, 2008)). The dominant deformation mechanism of this framework 
might involve a 'wine-rack'-like expansion in one direction that couples to a con- 
traction in a perpendicular direction (Baughman et al., 1998) — a mechanism that is 
reflected in the lattice parameter variation of the orthogonal (rectangular) unit cell 
shown in Fig. 1. In that of the related rhombic cell [Fig. 1(b)], there is no appar- 
ent change in the values of the lattice parameters a and b; instead the important 
variation occurs in the angle j3, which is less readily interpreted. By analogy, in real 
materials of low crystal symmetry the key science can be hidden within the variation 
of the entire ensemble of lattice parameters a, b, c, a, (3, 7, making analysis much less 
straightforward than would otherwise be hoped (or indeed is often assumed). 

The tensor algebra that resolves this issue is well established but is perhaps used 
more widely within some communities than others. The central concept is that for any 
system there exists a unique set of orthogonal axes (the principal axes) along which the 
material responds in a purely linear fashion — either expanding or contracting — with 
the response along other directions involving some non-zero shear component. It is 
the expansion or contraction along the principal axes that describes the fundamental 
mechanical response of the system. Predictably, crystal symmetry plays a crucial role 
in determining the orientation of these axes. For systems of orthorhombic symmetry or 
higher the principal and (conventional) crystallographic axes coincide (Nye, 1957). 2 

Hence the variation in lattice parameters I with temperature and pressure can be 

2 For uniaxial systems (i.e., with 3-, 4- or 6-fold rotational symmetry), two of the three principal axes 
are related by this rotational symmetry. These axes define a unique plane in which the mechanical 
response is invariant. Cubic systems are isotropic. 
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related directly to the bulk expansivities ae and compressibilities Kf 



.3 



K t = 



a e = + 




(2) 



(1) 



For monoclinic and triclinic crystal classes, however, the principal axes are not neces- 
sarily related to the unit cell axes 4 and must instead be determined as the eigenvectors 
of the full strain tensor. In particular, because the principal expansivities and com- 
pressibilities now depend in a complex fashion on the variation in unit cell lengths 
and angles, the values determined using Eqs. (1) and (2) no longer carry any physical 
significance. Returning to the example of Fig. 1, the discrepancy in lattice thermal 
expansivities arises because the principal axes coincide with the axes of the rectangular 
unit cell rather than the rhombic cell. 

Variable-temperature (or -pressure) unit cell data contain all the information nec- 
essary to identify the principal axes and hence the principal coefficients of thermal 
expansion (or compressibilities) — but for low-symmetry systems the process of calcu- 
lating these is a problem best suited to automation. A number of software packages 
are capable of carrying out this calculation: notably strain (Ohashi, 1982) and its 
successor winstrain, 5 deform (Filhol et al, 1987) and eIam (Marmier et al, 2010). 
Yet for much of the materials chemistry community, the many features offered by these 
packages can in some cases require more detailed experimental data than are available 
(e.g. elastic constant measurements), and/or provide a level of analysis beyond what 
is necessary to identify expansivity and compressibility anomalies. 

3 There is more than one convention widely in use for the symbols for compressibility and its inverse 
the bulk modulus. In key solid-state physics texts (e.g. Ashcroft & Mermin (1976) and Kittel (2005)) 
K is used to denote compressibility and B the bulk modulus. Confusingly, the primary reference 
for mechanical properties of crystals (Nye, 1957) uses K to denote the bulk modulus and (3 the 
compressibilities. Here, we follow the physics convention in order to remove any ambiguity regarding 
the lattice parameter /3. 

4 For monoclinic classes in the conventional setting, one principal axis lies parallel to b. 

5 The winstrain program is unpublished but available for download from www.rossangel.com. 
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We have developed a new web-based tool, PASCal (= principal axis strain calculatorl; 
available at http://PASCal.chem.ox.ac.uk) that accepts variable-temperature or variable- 
pressure lattice parameter data as input and simply returns as output the principal- 
axis expansivities/compressibilities and the orientation of the principal axes relative to 
the crystallographic axes. The pressure-dependence of the unit cell volume is also fit- 
ted to equations of state, from which the bulk modulus and its pressure derivative are 
determined. This output enables the user to identify quickly any anomalous mechan- 
ical properties (e.g. NTE, NLC, pressure-induced softening) and also to understand 
their geometric relationship to the structure itself. At all times, our particular focus 
has been on providing a simple tool useable by experts in elastic property calculations 
and non-experts alike. 

Our paper is arranged as follows: in Section 2, we describe the particular implemen- 
tation of strain tensor calculations used by the program and provide a summary of the 
most important aspects of the program input and output. What follows subsequently 
is a series of three case studies — focusing on one MOF, one molecular crystal, and 
one mineral — where we have used PASCal to re-analyse previously-published lattice 
parameter data. In all cases we are able to develop the relevant science as a result 
of considering the principal compressibilities or expansivities. In this way, we hope to 
demonstrate that routine application of PASCal to parametric crystallographic studies 
is not only feasible, but that the program is capable of adding substantive scientific 
value to the field. 

2. Computational Detail 

Theory 

The calculations at the heart of PASCal involve determination of the orthogonal 
strains. Our starting point is the transformation from crystallographic axes Aj to 
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orthogonal axes Ej (with notation as in Giacovazzo et al. (2011)). The orientation of 
the Ej is arbitrary: we make use of the Institute of Radio Engineers convention where 
E3 is parallel to the c crystallographic axis, Ei is parallel to a*, and E2 = E3 x Ei. 
The corresponding change-of-basis transformation is described by the square matrix 
M: 

E = M x A. (3) 

The strain e corresponding to a pair of initial and final conditions (i.e., temperature 
or pressure values) is given as the symmetric part of the product of M at the first 
condition with its inverse at the second condition: defining 

e = M final >< Minitial - I, (4) 

we have 

e= 1 -(e + e T ). (5) 

The eigenvalues and eigenvectors Xj of the matrix e are then the principal strains 
and the principal axes. The derivatives of these principal strains with respect to tem- 
perature and/or pressure give the principal coefficients of thermal expansion and/or 
compressibilities: 

Formally, the definitions we have given here correspond to infinitesimal Lagrangian 
strains — this is the standard approach applied in texts such as Nye (1957). We note 
however that infinitesimal strains are quantitatively exact only in the limit of small 
variations in lattice parameters. Consequently, slightly different definitions of strain 
can be more meaningful in other situations (e.g. finite strain for lattice discontinu- 
ities at a phase transition). There is scope within PASCal to calculate these alter- 
native strains; interested readers are referred to Giacovazzo et al. (2011), Schlenker 
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et al. (1978) and Zotov (1990) for further details regarding the choice of strain and 
the associated calculations. In practice the use of different strain definitions can give 
quantitatively significant differences in principal axis expansivity /compressibility val- 
ues, but is unlikely to affect the diagnosis of qualitative features, such as the existence 
or otherwise of NTE or NLC. Implementation of all strain calculations within PASCal is 
as a Fortran90 routine that makes use of the LAPACK and minpack libraries (Anderson 
et al, 1999; More et al, 1984). 

Program input 

We have tried to simplify the required input for PASCal as much as possible. Para- 
metric data are entered as a series of whitespace- or comma-separated values in 
sequential columns: temperature (K) or pressure (GPa) values; the uncertainty in 
these values; followed by the corresponding lattice parameters a, b, c, a, (3, 7. At this 
stage, estimated errors in cell parameters are neither required nor taken into account, 
as the errors in all subsequently calculated parameters are nearly always dominated 
by errors in temperature or pressure values and scatter in the data themselves. 6 In 
Fig. 2 we show an example set of input data for both expansivity and compressibility 

measurements of the high pressure phase Ag3[Co(CN)e]-II (Goodwin et al, 2008). 

6 The use of standard errors on lattice parameters to determine standard errors in derived param- 
eters has been the subject of a number of detailed investigations. First, it is well established that 
standard uncertainties are underestimated by both single-crystal and powder (Rietveld) refinement, 
with reported errors smaller by up to an order of magnitude than the real values determined by 
multiple measurements (Taylor & Kennard, 1986; Herbstein, 2000). Second, the avoidance of bias in 
derived orthogonal cells arising from angular uncertainty is a well-recognised problem for monoclinic 
and triclinic systems (Haestier, 2009). Third, the covariances necessary for accurate error propaga- 
tion are rarely quoted or included in deposited CIFs. Some of these issues can be overcome with 
careful treatment — e.g. by using repeat measurements to obtain experimental distributions of cell 
parameters — hence providing physically-meaningful estimates of standard uncertainties. A thorough 
investigation of one such approach is detailed in Angel (2000). But, in order to maintain a simple inter- 
face and in order not to give an impression of misleading precision, the error-propagation algorithms 
within PASCal have been implemented in such a way that is likely to provide a small overestimation 
of uncertainties. For specific instances demanding greater control over error treatment, the user is 
referred to more specialist code — e.g. strain (Ohashi, 1982). 
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Program output 

For variable-temperature data, PASCal calculates the principal thermal expansivities 
across the entire temperature range for which data are provided. No attempt is made to 
fit the expansivity data at any level of complexity beyond a simple linear dependence: 
£i(T) = £o(l + OLiT). More complex fitting approaches, such as using a Debye model 
for thermal expansion, are prone to instabilities in the absence of very high quality 
data that extend to low temperatures T <C 0d (David et al., 1993). We are of the 
view that complex fits are most likely to be of interest to specialist users, with simple 
linear coefficients of thermal expansion of greatest value to the broadest user base. 
In any case, the principal strains are given in tabular format for each temperature 
point, with the intention that more complex fits might easily be performed using any 
third-party fitting software of choice. 

The relationship between principal axes and the crystallographic axes is also described 
in the PASCal output, where we give the normalised components of the principal axes 
projected onto the crystallographic axes. The thermal expansivity tensor is repre- 
sented visually in the form of the expansivity indicatrix. The expansivity indica- 
trix is a smoothly-varying surface centred at the unit cell origin, where the distance 
between the surface and the origin in a given direction r is equal to the magnitude 
of a in that same direction (taking negative values as appropriate) (Belousov & Fila- 
tov, 2007; Durka et al, 2011; Krashnenko et al, 2000; Vershinin et al, 2005). Perhaps 
the more traditional method of visualising the expansivity or compressibility tensor 
would be via either the representation quadric or the strain ellipsoid (Nye, 1957). The 
particular advantage of the expansivity /compressibility indicatrix over these other 
representations — and the reason behind its inclusion in the PASCal output — lies in the 
ease with which crystallographic directions corresponding to mechanical behaviour of 
particular interest can be identified {e.g. a very large and positive, very large and 
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negative, or zero). We illustrate this point explicitly in Fig. 3, where we show together 
the expansivity indicatrix, the quadric and the strain ellipsoid for the one system 
Ag3[Co(CN)g]-II • For this material, comparison of the expansivity indicatrix with the 
crystal structure shows — in accord with the original study (Goodwin et al, 2008) — 
that NTE occurs along the layer-stacking axis and that the covalent Co-C-N-Ag- 
N-C-Co linkages are oriented along directions for which the expansivity essentially 
vanishes. The actual PASCal output for this same system is shown in Fig. 4. 

The output provided for variable-pressure data is more detailed. At the simplest 
level, median principal compressibilities, their uncertainties, and the corresponding 
principal axes are calculated according to Eq. 7 using a linear fit to the strains. 
While thermal expansion data are usually approximately linear across typical mea- 
surement ranges, a linear pressure-dependence is rarely observed for variable- pressure 
data (and thermodynamically forbidden over large pressure ranges). In our experience 
of relatively highly-compressible materials such as molecular frameworks, axial com- 
pressibilities are often much better characterised using an empirical fit of the form 
Zip) = ^Q + Xip — PcY (Goodwin et al, 2008; Cairns et al, 2012). Accordingly, PASCal 
also outputs fitted values of A and v for each principal axis and — because compressibil- 
ity now varies with pressure — the Ki (p) for each pressure point and their uncertainties 
are both plotted graphically and listed in tabular form [Fig. 5]. 

There is valuable thermodynamic information in the pressure-dependence of the unit 
cell volume, in the form of the bulk modulus B = (— Vdp/dV)T and its (dimensionless) 
pressure derivative B' = (dB/dpjx- The parameter B' can be an indicator of unusual 
mechanical behaviour; for example, a large B' is often indicative of a rapid stiffening 
characteristic of layered materials (Munn, 1972) and a negative B' can indicate the 
presence of a dynamic instability (Chapman & Chupas, 2007). Accordingly, unit cell 
volumes are fitted in PASCal using Birch-Murnaghan equations of state (Birch, 1947) of 
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various degrees of complexity, with the ultimate choice left to the user, perhaps guided 
by the quality of fits obtained or by construction of a so-called f-F plot (Angel, 2000). 
Defining the parameter 



where Vq is the zero-pressure unit cell volume, the second- and third-order Birch- 
Murnaghan fits correspond to the equations of state 

P(V) = 3 ^(r 1 7 -V 5 ) (9) 

and 

p(v) = ^ (V - v 5 ) 



1 + \{B> -4)^-1) 



(10) 



respectively. In Eq. (10), Bq can be interpreted as the value of B at zero pressure. For 
completeness, the option is provided to use finite Eulerian rather than infinitesimal 
Lagrangian strains, as the former is used in the derivation of the Birch-Murnaghan 
equation of state (Birch, 1947). We note, however, that this formalism was derived 
for isotropic crystals and need not describe equally well the behaviour of anisotropic 
systems. Consequently Birch-Murnaghan fits for anisotropic systems should be inter- 
preted with this caveat in mind. 

For high-pressure phases that are unstable below some critical pressure p c > 0, 
even the third-order Birch-Murnaghan equation of state does not capture the elastic 
behaviour fully, as it assumes that both V and B are continuous functions of p for all 
p > 0. In PASCal, we have incorporated the option to fit a modified version, which now 
takes into account the existence of a non-zero critical pressure p c (Sata et al, 2002): 

1 



Pc --(3B-5 Pc )(l-ri z ) 



p(V) = rf 



(") 



In fitting Eq. (11) the assumption is made that p c is less than or equal to the lowest 
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pressure value for which data are given. The user is given the option to include a fixed 
value of p c in their analysis. 

For the example of Ag3[Co(CN)e]-II discussed above and shown in Fig. 2, it is mean- 
ingful to fix the value of p c at the phase-I/II transition pressure of 0.19 GPa. A good fit 
to the experimental V(p) data is then obtained by PASCal, as shown in Fig. 5. The fun- 
damental parameters we extract in this fit show some discrepancy with those reported 
in Goodwin et al. (2008): here we obtain B = 14.2(14) GPa and B' = 10.9(13) (the 
previously-reported values are 11.8(7) GPa and 13.5(12), respectively). The origin of 
this discrepancy appears to be the extremely large parametric covariances within 
the third-order Birch-Murnaghan equation of state resulting in different convergence 
behaviour with different minimisation routines (the values reported in Goodwin et al. 
(2008) were determined using the solver routine within Microsoft Excel). A fit using 
the program EOS-FIT (Angel, 2001) (which automatically constrains p c = 0, and hence 
gives a slightly-different set of parameters yet again) returns covariance terms with 
magnitudes greater than 0.94, indicating that the system of equations is very poorly 
constrained. It is actually the case for many well-studied systems that there exists a 
variety of B, B' and Vb values that result in very similar fits to experimental compress- 
ibility data, and we reiterate here the caution of e.g. Jackson (1998) by emphasising 
the covariance amongst these parameters, especially for data with large scatter. 



3. Case Studies 

3.1. A metal-organic framework: Cu-SIP-3 

The first of our case studies concerns the flexible metal-organic framework Cu2(OH)(CsH307S)(H20)-2H20 
(abbreviated to Cu-SIP-3-3H20). This is a material of particular currency for its 
potential applications as a biomedical agent for controlled nitrous oxide (NO) release 
(Xiao et al, 2009; Hinks et al., 2010). Under ambient conditions its structure can 
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be considered as a stack of two-dimensional, covalently-connected 'sheets': each such 
sheet is assembled from what is technically a distorted square (4,4)-net of tetrameric 
CU4 clusters connected via 5-sulfoisophthalate linkers (although a description in terms 
of hexagonal close packing of the same clusters might be considered equally valid) 
[Fig. 6(a)]. Adjacent sheets interact only very weakly — via axial Cu 2+ . . .O3SR con- 
tacts (r(Cu. . .0) > 2.8 A) — such that if only strong covalent bonds are considered then 
the network of the solid as a whole has two-dimensional rather than three-dimensional 
connectivity. 

Variable-temperature (150-365 K) lattice parameter data for this ambient phase 
were reported by Allan et al. (2010) as part of a detailed study into a reversible, 
reconstructive phase transition that occurs on desolvation/resolvation. Because Cu- 
SIP-3-3H20 crystallises in the monoclinic space group P2i/n, this is an example of a 
system where not all of the coefficients of thermal expansion determined from lattice 
parameters correspond to principal expansivities. Nonetheless there are hints even 
in the lattice parameter measurements that the system exhibits anomalous thermal 
expansion behaviour: the b unit cell parameter decreases on increasing temperature, 
and there is a negligible change in unit cell volume over the temperature range studied. 

We have used the lattice parameter data of Allan et al. (2010) as input for PASCal 
in order to determine the fundamental thermal expansion properties of this interest- 
ing material; our results are summarised in Table 1. Perhaps the most unexpected 
finding here is that two of the three principal expansivities are in fact negative. So 
not only does Cu-SIP-3-3H20 become one of the few known materials to show an 
area-NTE effect, but this effect is very large indeed. In fact the area coefficient of 
thermal expansion a a = ol\ + «2 = — 76MK -1 is essentially an order of magnitude 
more negative than corresponding values in previously-studied area-NTE materials: 
e.g. the Ni(CN) 2 family (a A = -11 to -16 MET 1 ) (Hibble et al, 2007; Hibble et 
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al., 2011), graphite (a^ = -2.4MK" 1 ) (Bailey and Yates, 1970) and low-temperature 
arsenic (a A = -4.0MK" 1 ) (Munn, 1972). 

Clues as to the microscopic driving force responsible for this unusual behaviour can 
be found in the orientation of the principal axes relative to the crystal structure. The 
set of directions in which NTE is observed lies perpendicular to the X3 principal axis 
(i.e., the plane containing the xi and X2 axes). Hence the orientation of X3 must carry 
some special significance. We find that X3 is approximately parallel to the [101] crystal 
axis, which is in fact the layer-stacking direction (i.e. normal to the layer shown in 
Fig. 6(a)). Hence, on heating, the framework Cu-SIP-3 expands in a direction parallel 
to the stacking axis, but contracts in every perpendicular direction. 

The general propensity for layered materials to exhibit area-NTE was in fact pre- 
dicted sixty years ago (Lifshitz, 1952): the lowest-energy phonon branches correspond 
to 'rippling' of the layers, a displacement pattern that causes the interlayer separation 
to increase and the effective layer area to decrease [Fig. 6b] . The actual magnitude of 
NTE response observed then depends on the elastic stiffness of the sheets (i.e., how 
easily they might be corrugated). We would suggest that the use of flexible organic 
linkers is responsible for the large NTE effect we find here by virtue of the low energy 
cost associated with flexing of the layers. Finally we comment that in the light of 
an emerging correspondence between negative thermal expansion and negative com- 
pressibility in framework materials (Goodwin et al, 2008; Fortes et al, 2011; Cairns 
et al, 2012), one might expect Cu-SIP-3 to exhibit the elusive property of negative 
area compressibility under hydrostatic pressure (Baughman et al, 1998). 

3.2. A molecular crystal: HMX 

Anisotropy in thermal expansion can be of considerable importance in materials 
applications, particularly when the avoidance of cracking is important. One such appli- 
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cation is the development of high-performance explosive materials, where cracking can 
increase the propensity for accidental detonation — i.e. the sensitivity of the explosive 
(Bennett et al, 1998; McGrane et al, 2009; Ramaswamy et al, 1996). Like many 
molecular crystals, explosives often crystallise in low-symmetry crystal systems. Con- 
sequently in order to determine the degree of anisotropy in their mechanical response 
it is rarely sufficient to consider the temperature-dependence of the unit cell lengths 
alone. Indeed since the principal coefficients of thermal expansion correspond to the 
extrema of mechanical response, the use of unit cell length variation to characterise 
mechanical anisotropy will always underestimate its maximum extent. 

One of the most important and widely- used explosives is cyclotetramethylene-tetranitramine 
(HMX) (Agrawal, 2010). HMX is known to crystallise in any one of four polymorphs 
(a, (3, 7, 5), of which the /3 phase (space group P2\/n) is the most important because 
it is the only stable polymorph under ambient conditions [Fig. 7(a)] (Cady et al., 1963; 
Choi and Boutin, 1970). A recent thermal expansion study of /3-HMX, aimed at charac- 
terising its mechanical anisotropy, revealed a large difference in the lattice coefficients 
of thermal expansion: a a ,aj„a c = 33,99 and 0.43 MK -1 , respectively (Deschamps et 
al, 2011). 

Using the same variable-temperature lattice parameter data of Deschamps et al. 
(2011) as input for PASCal, we have been able to determine the principal coefficients 
of thermal expansion; the PASCal output is summarised in Table 2. We find a sub- 
tle but important difference between the principal and lattice coefficients of thermal 
expansion: namely that there exists NTE along one of the principal axes (X3). This 
is important because the degree of anisotropy in some directions is now substantially 
larger than anticipated from the lattice parameters alone; for example, the lattice 
expansivities underestimate anisotropy in the (010) plane by ca 40%. 

At a molecular level, the driving force for the NTE effect we observe appears to be a 
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flattening of the eight-membered ring of each HMX molecule at higher temperatures. 
There are two different ring orientations within the /3-HMX crystal structure, which 
are related by the space group symmetry. For each orientation the four ring N atoms 
are coplanar, so a vector normal to this plane characterises the ring orientation. We find 
that an average of the two orientation vectors coincides with the NTE principal axis 
X3 to within 1% [Fig. 7(b)]. A simple measure of the degree of ring puckering is given 
by the mean-squared displacement of the four C atoms away from the N-atom plane 
defined above. Based on the atomic coordinates given by Deschamps et al. (2011), this 
quantity appears to decrease with temperature — i.e. the HMX rings become flatter at 
higher temperatures — providing a feasible explanation for the reduction in crystallite 
size in a direction perpendicular to these C4N4 rings. 

3.3. A mineral: malayaite 

Our final case study concerns the mineral malayaite (CaSnOSiC^), which is a 
calcium tin silicate closely related to titanite (CaTiOSi04) that — when doped with 
chromium — finds application as a pink glaze for ceramics (Ingham & Bradford, 1961; 
Higgins & Ribbe, 1976). From a fundamental science viewpoint, the comparative crit- 
ical behaviour of the malayaite-titanite family is of general interest in helping under- 
stand the interplay between mechanical and electronic properties in framework sili- 
cates (Salje et al., 1993; Zhang et al, 1999). In titanite itself, second-order Jahn Teller 
distortions of the d° Ti 4+ ion lead to local polarisations of the TiC>6 octahedra (such 
as is observed in key ferroelectrics such as BaTiC>6 (Shirane et al, 1967)); this effect is 
reduced on substitution with Sn 4+ . An interest in understanding the effects of pressure 
and temperature on the interplay between these off-centre distortions and structural 
transitions has led to a number of variable-temperature and variable-pressure studies 
on both systems (Kek et al., 1997; Kunz et al, 2000). 
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The structure of malayaite consists of chains of vertex-sharing SnOg octahedra. 
These chains are connected via SiC>4 tetrahedra to form a three-dimensional anionic 
network structure. Calcium ions occupy a set of voids within this framework such that 
each Ca 2+ is 'coordinated' by six or seven O atoms (the Ca-0 separations being suffi- 
ciently large that identification of the coordination sphere is somewhat arbitrary) . At a 
pressure of 4.95(1) GPa, malayaite transforms from its ambient monoclinic A2/a phase 
into a triclinic PI phase (Rath et al., 2003). The distortion mechanism associated with 
this transition and also the structural changes observed on further compression within 
the triclinic phase offer insight into the key dynamical properties of malayaite. So, for 
example, the high-pressure crystallographic study of Rath et al. (2003) documents, 
amongst other things, the effective compressibilities of SiC-4 tetrahedra (not very com- 
pressible), SnC-6 octahedra (more compressible) and C&O x coordination polyhedra 
(very compressible); the implication being that vibrational modes resulting in large 
changes in Ca-0 distances play a significant role in the dynamics of ambient-phase 
malayaite. 

Making use of the variable-pressure lattice parameter data reported by Rath et al. 
(2003) for the high-pressure triclinic phase of malayaite, we have been able to extend 
their analysis to the determination of orthogonal strains. Using first a linear fit to 
these data, PASCal obtains the principal compressibilities and principal axes given in 
Table 3. Of particular note is the NLC behaviour along X3 — i.e., approximately par- 
allel to the [441] crystallographic axis. We find that the compression mechanism is 
dominated by concerted rigid-body rotations and translations of the Sn06 and Si04 
polyhedra that cause a contraction within the plane perpendicular to the [441] axis. 
The polyhedral 'rocking' associated with this compression translates to an expansion 
in a perpendicular direction — the NLC axis. In this process the CaO^ coordination 
environment undergoes substantive changes in geometry, such that one might sug- 
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gest the coordination flexibility of Ca 2+ allows a particularly anisotropic mechanical 
response to network compression. 

For completeness, we show in Fig. 8 the non-linear parameterised fits to the princi- 
pal strains obtained using PASCal and the derived compressibilities. What is apparent 
from these plots is that the value of K2 increases with pressure — even in this post- 
transition phase — which could suggest the existence of an incipient dynamic instability 
at yet higher pressures. This conclusion is reflected also in the thermodynamic quan- 
tities obtained from a third-order Birch- Murnaghan fit to the V(p) data: here we find 
Bo = 170(40) GPa and B' = —3(5). Given the magnitude of the uncertainties involved, 
we cannot be sure that the value of B' is really negative, but if it were then the impli- 
cation is that the material is becoming mechanically softer on compression — a feature 
also associated with dynamic instabilities (Chapman & Chupas, 2007; McConnell et 
al, 2000). To some extent this was discussed by Rath et al. (2003), but caution led the 
authors to constrain their fit to a second-order Birch-Murnaghan equation of state. If 
we perform a similar fit, PASCal obtains a value of B = 119(2) GPa that is in close 
agreement with the value obtained by Rath et al. (2003) (B = 118.3(7) GPa). Intrigu- 
ingly, Rath (2002) reports indications of a second high(er)-pressure phase transition 
in titanite at ca 10 GPa, suggesting that the dynamical anomalies we observe here 
for malayaite might signal the existence of a similar incipient transition. This puta- 
tive phase transition could provide an explanation for the observed NLC, as has been 
observed in other oxide frameworks (Angel & Bismayer, 1999; Angel et al, 2004). 



4. Conclusion 

Our aim has been to develop a simple tool that enables rapid and straightforward 
calculation of the principal coefficients of thermal expansion and compressibilities from 
variable-p/T lattice parameter data. On the one hand, this will allow users to screen 
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easily for anomalous mechanical responses — e.g., NTE, NLC and pressure-induced 
softening. On the other hand, the reporting of principal axes means that any such 
effects can be related back to the underlying crystal structure. While we are aware 
that PASCal replicates some functionalities of existing codes, it is our hope that the 
straightforward input format, web accessibility, and ease of output interpretation will 
mean that a greater scientific audience can make use of orthogonal strain calculations 
than is currently the case. 

The three case studies discussed above illustrate some of the ways in which scientific 
value might be added through routine lattice parameter analysis with PASCal. In the 
case of the metal-organic framework Cu-SIP-3, we have shown it exhibits the strongest 
area-NTE effect yet reported, and have been able to relate this effect back to low- 
energy acoustic modes of its layer-like structure. Analysis of thermal expansion data for 
the explosive material HMX reveals a much stronger mechanical anisotropy than had 
otherwise been determined: a result that may have implications for strain propagation 
and cracking of polycrystalline HMX preparations. Finally, we have shown that high- 
pressure malayaite exhibits two thermodynamic anomalies: it expands in one direction 
on compression, and it also becomes mechanically softer on increasing pressure. 

These three systems are unlikely to be isolated examples. Indeed, we anticipate that 
with a more widespread application of principal axis strain calculations, phenomena 
such as NTE and NLC will likely be observed with increasing frequency. Irrespective of 
the prevalence of such effects, the mechanisms responsible will almost always provide 
important insight into the dominant dynamical properties of materials. So it is our 
particular hope that the ready availability of programs such as PASCal might help 
develop our understanding of dynamics in classes of solid materials for which strain 
calculations have not yet played a central role. 
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Table 1. Principal coefficients of thermal expansion and corresponding principal axes 



determined for Cu-SIP-3-3H 2 0. 







Component of Xj along the 




Principal 




crystallographic axes 




axis, i 


(MK- 1 ) 


a be 


Approximate axis 


1 


-55(7) 


-10 


[010] 


2 


-21(7) 


-0.7932 0.6090 


[101] 


3 


+60(9) 


0.8260 0.5637 


[101] 



Table 2. Principal coefficients of thermal expansion and corresponding principal axes 

determined for /3-HMX. 







Component of Xj along the 




Principal 




crystallographic axes 




axis, i 


(MK" 1 ) 


a be 


Approximate axis 


1 


+101(3) 


-10 


[010] 


2 


+37(3) 


-0.5570 0.8305 


[203] 


3 


-15(3) 


-0.8419 -0.5397 


[302] 



Table 3. Principal compressibilities and corresponding principal axes determined for the 
high-pressure (PI) phase of malayaite. 







Component of Xj along the 




Principal 


Ki 


crystallographic axes 




axis, i 


(TPa- 1 ) 


a b c 


Approximate axis 


1 


+7.6(4) 


-0.6482 -0.5429 0.5339 


[111] 


2 


+2.28(10) 


-0.5846 -0.0820 -0.8072 


[203] 


3 


-3.1(3) 


-0.6270 0.7591 0.1747 


[441] 
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Fig. 1. An illustrative example of the effect of unit cell choice on lattice parameter 
expansivity determination: (a) A hypothetical two-dimensional molecular frame- 
work, with both a rhombic and an orthogonal unit cell indicated; (b) a plot of the 
relative change in lattice parameter determined for each of these two cells for the 
concerted lattice deformation indicated by arrows in (a). The rhombic cell shows no 
variation in unit cell length whereas the changes in orthogonal cell lengths reflect 
the underlying bulk material response. 
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A web tool for Principal Axis Strain Calculations 

Quick Guide . Technical information and caveats . Mote information available or 



Phase transition: 

Use a critical pressure?: 

Critical pressure: 

Advanced Options: 
Finite strain? 
Eulerian strain? 



(a) 



Developed by Matthew Cliffe and Andrew Goodwin at the University of Oxford. 
Any queries or comments to matthew[dot]cliffe[at]chem[dot]ox[dot]ac[dot]uk. 
Please cite M.J. Cliffe and A.L Goodwin. arxiv.orq:1204.3007 . 




A web tool for Principal Axis Strain Calculations 
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Data type: 

Temperature (K) Pressure (GPa) 

Phase transition: 

Use a critical pressure?: a 

Critical pressure: 1 0.19 



Advanced Options: 
Finite strain? 
Eulerian strain? 

Data: 

TorP o(TorP) a b c a p y 




(b) 



Developed by Matthew Cliffe and Andrew Goodwin at the University of Oxford. 
Any queries or comments to matthew[dot]cliffe[at]chem[dot]ox[dot]ac[dot]uk. 
Please cite M.J. Cliffe and A.L. Goodwin. arxiv.orq:1 204.3007 . 



Fig. 2. Sample PASCal input for (a) variable-temperature and (b) variable- pressure 
lattice parameter data measured for Ag3[Co(CN)e]-II (Goodwin et al, 2008). 
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(a) (b) (c) 



Fig. 3. A comparison of three methods for representing an expansivity tensor, using the 
variable-temperature Ag3[Co(CN)e]-II data of Goodwin et al. (2008). (a) The ther- 
mal expansivity indicatrix, with PTE shown in red and NTE in blue; (b) the thermal 
expansivity quadric; (c) the strain ellipsoid (shown here for the maximum strain 
obtained in order to emphasise the difficulty associated with perceiving aspheric- 
ity). The directions and relative magnitudes of large PTE and NTE are obvious 
in the expansivity indicatrix representation, but are obscured in the other, more 
traditional, representation surfaces. 
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Fig. 4. (a) The key PASCal output obtained for the input data shown in Fig. 2(a), 
showing the principal components of the expansivity tensor and the expansivity 
indicatrix representation as discussed in the text, (b) A representation of the crys- 
tal structure of Ag3[Co(CN)g]-II shown with axes in the same orientation as the 
expansivity indicatrix in (a). The layer of Ag + cations (shown as large spheres) cor- 
responds to the large positive thermal expansion component of the expansivity indi- 
catrix, which is shown in red; the perpendicular (stacking) axis of Ag3[Co(CN)e]-II 
then corresponds to the large NTE component of the expansivity indicatrix (shown 
in blue). The Co-C-N-Ag-N-C-Co linkages, shown in ball-and-stick representa- 
tion, are oriented in the same set of directions for which the expansivity indicatrix 
has a near-zero value. 
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Fig. 5. Calculated fits to variable-pressure Ag3[Co(CN)e]-II lattice parameter data 
obtained using PASCal: (a) relative changes in length for the three principal axes, 
(b) the corresponding principal compressibilities, and (c) the modified third-order 
Birch-Murnaghan fit to the unit cell volumes. 
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Fig. 6. Negative thermal expansion in Cu-SIP-3-3H20. (a) A single covalently-bonded 
layer in the material, viewed down the [101] axis (the single positive thermal expan- 
sion axis), (b) The general mechanism responsible for area-NTE in layered mate- 
rials: the lowest energy deformations correspond to 'rippling' of the layers, which 
results in an increased interlayer separation but a smaller effective layer area. 
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Fig. 7. (a) A representation of the molecular geometry of HMX in its /3-phase. (b) 
The thermal expansivity indicatrix viewed down the axis of NTE (i.e. for the xi, X2 
plane). The negative- valued lobe of the expansivity indicatrix is small and projects 
from the centre of the expansivity indicatrix towards the viewer, (c) A represen- 
tation of the crystal structure of /3-HMX viewed in the same orientation as for 
the expansivity indicatrix in (b). A flattening of the central eight-membered C4N4 
rings at higher temperatures causes an expansion in the xi and X2 directions but a 
contraction in the X3 direction (i.e. along the viewing direction). 
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Fig. 8. Compressibility of high-pressure malayaite. (a) Relative changes in the lengths 
of the principal axes (open circles, solid lines) and PI unit cell axes (crosses and 
dashed lines) under compression. The variation in unit cell axes is substantially 
more moderate than the variation in principal axes, (b) The corresponding principal 
linear compressibilities. The increase in the magnitude of if 2 observed at increasing 
pressure reflects a dynamical instability. 
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Synopsis 

The web-based program PASCal (http://pascal.chem.ox.ac.uk) is introduced as a tool for deter- 
mining principal coefficients of thermal expansion and compressibilities using variable-pressure 
and variable-temperature lattice parameter data. The additional scientific insight provided by 
this analysis is illustrated with reference to three case studies: the metal-organic framework 
Cu-SIP-3, the molecular crystal /3-HMX and the mineral malayaite. 
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